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Abstract 


equation for calculating the pressure drop gas flow accompanied 
large temperature change, derived from the differential form the Bernoulli 
theorem combined with the heat transfer rate equation, was examined experi- 
mentally. Heating runs were made the Reynolds number range 6000 
15,000, for air flowing through single horizontal in. standard steel pipe with 
steam jacket. pressure drops calculated from the proposed equation were 
found agree very closely with the observed experimental values. The method 
evaluating pressure drop for turbulent nonisothermal flow assuming 
isothermal conditions temperature midway between the surface temperature 
and the mean bulk temperature the fluid, first developed McAdams for 
certain oils, was also tested for number selected runs. This method was 
found accurate the previous one, and have the advantage being 
less cumbersome application. 


Introduction 


The problem fluid friction for isothermal flow round pipes has been 
very thoroughly investigated compared with that for nonisothermal flow. 
The friction factor defined (16) the familiar differential form 
Bernoulli’s theorem: 

udu 


where the symbols have the meanings listed the nomenclature. For stream- 
line viscous flow, the law established Hagen (4) 1839 and Poiseuille 
(11) 1842 has been many times confirmed the calming length were 
least diameters terms the present symbols 

For turbulent flow, the accepted friction factor plots for smooth tubes and 
clean commercial iron and steel pipes, respectively, are those Drew, Koo, 
and McAdams (3), the results all-embracing literature survey. the 
case corroded tuberculated commercial pipes, however, especially those 
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smaller diameter the present investigation, values are considerably 
greater than those read from the aforementioned plots (2, 7), and roughness 
factor must applied. 


Recorded research nonisothermal fluid friction has far been restricted 
liquids only, particularly oils. Available data (1, 18) all indicate 
that for given value the Reynolds number based viscosity evaluated 
the mean bulk temperature the oil, the friction factor during heating 
less than isothermal flow, while, for cooling, the friction factor greater, 
especially when the flow the viscous region. These results have been 
reconciled with the isothermal friction factor plots evaluating the viscosity 


temperature 
rad 


t)/4 


for streamline flow, and 


of 


for turbulent flow (10). Another successful correlation, based the same 
data well other data (14), involves multiplying the friction factor for 
nonisothermal flow correction where 


the turbulent region study the cooling water brass 
tube (12), was found that nonisothermal friction factor data fall line 
with isothermal data when the viscosity term the Reynolds number 
evaluated the mean wall temperature. complex triple integral expression 
(8, pp. 571-577) relating isothermal and nonisothermal pressure for 
viscous flow has been verified for water and oil (9). 


1.1 
1.0 


for viscous flow, and 


the case nonisothermal flow gases, since the variation viscosity 
with temperature small compared with that liquids, the problem 
finding the correct temperature which evaluate order use iso- 
thermal friction data minor consequence (6). But the compressible 
property gases introduces the new problem how evaluate correctly 
the changing specific volume, which occurs explicitly the pressure head term 
Bernoulli’s equation, and implicitly the velocity head and Fanning 
friction terms. Where the specific volume does not vary greatly, the 
Bernoulli differential equation may integrated with little error using 
average value (17). However, this not applicable cases where 
the variation large. 


method integrating the differential form the Bernoulli equation 
assuming that the pressure change small compared with the temperature 
change was developed Holden and Gauvin McGill University 
1943, using simple thermodynamic relations. The present investigation 
was undertaken for the purpose testing this formula within the lower 
turbulent range for the heating air. 
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Formula Derivation 


The energy balance for the flow gases has been developed follows (17). 
The differential form Bernoulli’s equation given 


udu 
vdp 


(1) 
for horizontal conduit. 


Converting from local velocity average velocity since and 
for steady flow, 


Multiplying through 
G*pdv 


Substituting the ideal gas law relation 


(2) 


For isothermal flow, this equation can integrated directly the familiar 


form 


For nonisothermal flow, Equation (2) rearranged: 


pdp 


Since the Reynolds number does not vary greatly, may considered 
constant, and average value used. Holden and Gauvin considered that 
the change pressure relatively small compared with the change 
temperature, then may substituted for and this leads 


expression that could integrated: 
pap Deg 


heat balance differential length dL, 


4hdL 
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Therefore 
Substituting the energy balance, 
2 2 


For rough pipe, must modified the roughness factor and the final 
expression becomes 


a L 


Description Apparatus 


diagrammatic sketch the apparatus used this investigation shown 
Fig. steady flow air was provided means rotary compressor, 
belt-driven hp. motor. The air rate was controlled means 
by-pass valve the atmosphere and was measured orifice. The 
orifice drop was read water manometer and the orifice upstream pressure 
mercury manometer. The steam pressure was regulated means 
globe valve and measured Bourdon gauge. calming length 
diameters preceded the heating section 5.0 ft.), which consisted 
in. standard steel pipe 0.0303 ft.) inside in. standard steel steam 
jacket, with (‘‘Aircell’’) lagging in. thick) the outside. 
Thermometers were placed the inlet and outlet the heating section, and 
the pressure drop across this section was read mercury manometer. 


q 
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STEAM 


PIPE 


MERCURY 
MERCURY AIR. FROM 
COMPRESSOR 
5.0 FT. 


Experimental Procedure 


order determine the roughness factor and its variation with mass 
velocity, set isothermal trials were made initially throughout the 
entire velocity range the apparatus. 

Four sets heating runs each were then made steam pressures 
10, 15, and p.s.i. gauge, respectively. run was carried for 
min. after ensuring that thermal equilibrium had been reached, and read- 
ings orifice drop and pressure, inlet, and outlet temperatures, and pressure 
drop, were recorded five-minute intervals. 

The mean values the above readings for each run were used Equation 
(5) calculate the pressure drop. The temperature corresponding 
saturated steam the recorded steam pressure was taken the inside 
surface temperature (¢,) the conduit, calculations showed that the air 
film was controlling. The calculated and observed pressure drops for the 
runs are given Table 


TABLE 


COMPARISON OF CALCULATED AND OBSERVED PRESSURE DROPS 


ts, G, AP cate. APors. 
687.1 2.92 1.44 


STEAM 
STEAM JACKET ORIFICE 
|STEAM TRAP 
WATER 
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TABLE I—Concluded 


COMPARISON CALCULATED AND OBSERVED PRESSURE 


te, G, APobs. ’ 

4.46 3.13 3.15 —0.63 

5.10 4.01 —1.49 

5.20 4.14 4.12 +0.48 

5.71 4.99 5.04 —0.99 

5.94 5.26 5.31 

6.14 5.63 5.66 

6.26 5.94 5.91 +0.51 

6.45 6.20 6.08 +1.97 

709.8 2.96 1.43 —5.30 

3.80 2.30 2.44 

3.95 2.45 2.60 

4.45 3.12 —2.80 

4.91 3.74 3.87 —3.36 

5.36 4.45 4.57 —2.62 

5.75 5.05 5.16 —2.13 

6.06 5.56 5.61 —0.89 

6.28 5.78 6.00 —3.67 

6.51 6.37 6.37 0.00 

718.8 2.86 1.36 

3.51 1.99 2.02 

4.51 3.22 3.24 —0.62 

4.91 3.76 3.82 

5.42 4.58 4.62 —0.87 

5.53 4.66 4.79 —2.71 

5.91 5.43 —1.84 

6.08 5.70 —1.05 

6.43 6.25 6.29 —0.64 

Discussion 


The results indicate substantial agreement between calculated and observed 
pressure drops for the entire range velocity and steam pressure used the 
investigation. Although the deviation range —6% within the 
possible experimental error, with the percentage deviations 
numerically decreasing with increase manometer readings would 
expected, the preponderance negative over positive deviations may not 
merely coincidence, and there were apparent apparatus defects 
account for it. The probable explanation lies the fact that the calculations 
for these runs were made using friction factor based the mean air 
temperature, rather than mean the friction factors the inlet and outlet 
the heating section. The latter procedure would have resulted slightly 
higher value each case, hence slightly higher calculated pressure drop, 
and therefore slightly higher (algebraically, not numerically) deviation. 
This contention was actually verified for few sample runs. There would thus 
have been more symmetrical distribution positive and negative errors. 
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should noted that, from the derivation Equation (5), this relation 
applicable only when the pressure change relatively small compared 
with the temperature change, and this disparity often the case industrial 
practice, e.g., air heaters and coolers, gas preheaters, etc. Although the 
present experiment has verified this equation only for the heating air 
the lower turbulent region (Re 6000 15,000), there appears nothing 
the derivation that would preclude its applicability the heating cooling 
any reasonably gas either the viscous flow higher turbulent 
ranges, the above condition assumed hold. However, not valid 
make such generalization without further experimentation. 

was deemed interest calculate the pressure drop for few the 
runs this investigation means the familiar relation (Equation (3) 
for isothermal pressure drop for gas. Since the present case the pressure 
change relatively small, and must modified the roughness factor 


the above equation reduces 


this calculation were those the lowest and highest mass velocities for each 
the four steam pressures. The results are tabulated Table 


TABLE 


EVALUATION DATA USING “ISOTHERMAL” EQUATION 


687.1 2.92 646.0 1.41 1.44 —2.08 —5.55 
6.40 649.5 6.10 5.98 +2.01 +0.33 
2.74 1.27 1.30 —2.31 —5.38 
6.45 656.8 6.30 6.08 +3.62 +1.97 
709.8 2.96 661.0 1.48 1.51 
6.51 664.3 6.49 6.37 +1.88 0.00 
718.8 2.86 668.8 1.40 1.36 +2.86 
6.43 669.5 6.41 6.29 —0.64 


The calculated results are therefore generally accurate those obtained 
using the much more complex relation Equation semiempirical 
method seemingly would more useful practice, least for the lower 
turbulent range. 

The fact that the deviations using Equation (3) are consistently greater 
algebraically than the corresponding deviations using Equation (5), and that 
there more symmetrical distribution positive and negative errors 
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the case the former, supports the contention that the uneven distribution 
deviations using the latter formula not due divergence the actual 
apparatus from ideal system. 


Conclusions 


Equation (5), developed Holden and Gauvin, successfully predicts, 
with reasonable accuracy, the pressure drop for the flow gas horizontal 
pipe accompanied large temperature change, within the lower turbulent 
region (Re 6000 15,000). 


The method evaluating data for turbulent nonisothermal runs assum- 
ing that these runs occur isothermally temperature midway between the 
surface temperature and the mean bulk temperature the fluid, method first 
developed McAdams for certain oils, accurate the method 
Holden and Gauvin predicting nonisothermal pressure drop gas 
flowing horizontal pipe, and has the advantage for industrial design 
purposes being much less cumbersome. 


Definition Symbols 


inside diameter conduit, ft., 

friction factor Fanning equation, dimensionless, 

acceleration due gravity, 32.2 

gas film coefficient heat transfer obtained from heat balance 


the heating section, B.t.u./(sec.) abs.), 
heating (or cooling) length conduit, ft., 
reciprocal molecular weight gas, 


calculated pressure drop, cm. Hg. 
(144) (14.7) 
observed pressure drop, cm. Hg., 


steam gauge pressure, 


; 

\ 

3 

og 

= 

» 
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ideal gas law constant, 1543 abs.), 
Reynolds number, dimensionless, DG/0.000672 
gas temperature, abs.), 

inlet temperature gas, abs.), 

outlet temperature gas, abs.), 

mean bulk temperature fluid, abs.), 

surface temperature, abs.), 

(t, t)/4 for below 2100, abs.), 

local velocity fluid, ft./sec., 

average velocity fluid, ft./sec., 

specific volume fluid, 


roughness factor, dimensionless, 


AP iheor. 


height conduit above arbitrary datum, ft., 
viscosity fluid mean bulk temperature, centipoises, 


viscosity fluid surface temperature, centipoises. 
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THE DETERMINATION THE EFFICIENCY HYDRAULIC 


JoHN KATZMAN 


Abstract 


Western Electric bead type thermistors were used determine the difference 
temperature the water between the headrace and tailrace hydraulic 
turbine different loads the turbine. From this temperature difference the 
maximum efficiency the turbine was determined 92% and the unit flow 
water cubic feet per second was calculated. 


Introduction 


The calorimetric method for the determination the efficiency 
hydraulic turbine was first developed Europe (1), but only small 


extent. The temperature water increased approximately 778 


when the water falls through height feet. When this water passes 
through hydraulic turbine the lower level, the amount energy taken 
from the water will show decrease the temperature the water; this 
decrease will depend upon the load the turbine. The difference the 
temperature the water the headrace and the tailrace was obtained 
means thermistors. 


This problem was brought the attention the author the Shawinigan 
Water and Power Company. The measurements were made the Shawinigan 
Water and Power Company powerhouse Quebec, Canada. 


Experimental 


Western Electric bead type thermistors whose resistance the ice 
point differed only 0.6% were chosen from group eight. These 
thermistors were chosen because the measurements were made when the 
temperature the water was The thermistors were fixed means 
paraffin wax brass tubes with small copper tips, Fig. small amount 
mercury was placed the bottom the tip for better thermal contact. 
The tubes were fitted with in. brass unions with which the unit was attached 
the supporting pipes that held them place the water. 

Each thermistor was connected one arm Wheatstone bridge, Fig. 
The other two arms were respectively, 10,000 ohm Leeds Northrup 
“unmounted and similar one 9500 ohms connected Leeds 
Northrup decade box 999.9 ohms. galvanometer (sensitivity, 0.26 
per mm: one meter) was used null point instrument. Six volts was 
impressed the bridge. The galvanometer was mounted shock proof 
support. 
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calibrate the bridge the thermistors were placed vessel two 
compartments, Fig. Each compartment measured ft. and was 
filled with 200 crushed ice that was intimately mixed with water. Both 

PLATINUM 


RESISTANCE 
THERMOMETER 


THERMISTOR — 


Fic. baths. 


thermistors were placed one compartment and bridge balance obtained. 
The thermistors were then placed separately the ice baths and platinum 
resistance thermometer was fixed beside one the thermistors. Small 
quantities sodium chloride were added the mixture which 
the resistance thermometer was placed. was thus possible obtain 
bath whose temperature could varied small amounts. For each change 
temperature the one bath the bridge was balanced. The change 
resistance necessary balance the bridge was plotted against the change 
temperature, Fig. 


Temperature difference, °F. 


1 


3.0 4.0 
Ohms 


Fic. curve the Wheatstone bridge arrangement. 


= HP D.C. 
STIRRER 
1.0 2.0 5.0 6.0 7.0 
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Tests 


The No. unit the nine units the powerhouse Quebec, 
was chosen for test. All tests were made Sundays because only then 
could the load the turbine varied will and the unit could put 
independent hand control. 


Both thermistors were lowered into the water the headrace the 
emergency gate opening and balance the bridge was obtained; this value 
was used the zero. Then one thermistor was placed position the 
headrace and the other the tailrace. The thermistor the headrace was 
moved and down for distance ft. determine whether the tempera- 
ture the water was uniform with depth. change temperature with 
depth was found. The same held for the tailrace. Temperature differences 
were measured for various gate openings, and the same time readings 
the turbine output and head water were obtained from the control room. 


Results 

The results obtained are shown Table Examination the divided 
differences the temperature difference indicated that, apart from random 
fluctuations, the relation the temperature difference the turbine output 
over the range studied could adequately represented parabolic equation. 


TABLE 
Turbine Gate Operating Gate Operating 
difference, output, head, ft. difference, output, ft. 

0.104 1.0 84.2 0.030 11.8 55.0 84.2 
0.087 85.4 0.017 84.1 
0.082 25.0 85.4 0.015 55.0 85.3 
0.087 2.0 25.0 84.3 0.012 12.6 85.1 
84.5 0.017 12.8 60.0 82.9 
0.052 3.4 85.3 0.012 13.4 69.0 85.4 
0.059 84.3 0.014 14.0 85.2 
4.2 30.0 85.5 0.007 15.0 65.0 85.4 
0.054 4.8 32.5 84.4 0.008 15.8 67.5 84.0 
0.042 0.007 15.8 67.5 85.1 
0.049 5.8 84.2 0.005 15.8 70.0 85.1 
0.042 6.0 35.0 85.4 16.8 72.5 
0.042 6.6 84.2 0.005 75.0 85.4 
0.025 6.8 85.2 17.6 85.1 
7.6 40.0 84.2 0.010 17.6 80.0 
0.027 7.8 40.0 85.4 0.007 17.8 80.0 85.3 
0.032 8.0 84.1 0.010 17.8 82.5 83.9 
0.020 8.6 42.5 85.1 0.010 18.0 82.5 85.1 
0.030 9.0 45.0 84.2 0.012 18.0 85.0 83.9 
0.018 9.4 45.0 85.4 0.010 18.2 87.5 83.9 
0.023 9.8 47.5 84.1 0.010 18.2 91.25 83.9 
0.015 9.8 47.5 85.1 0.009 18.3 85.0 85.3 
0.027 10.6 50.0 84.2 0.012 18.3 87.5 85.0 
0.022 10.8 50.0 85.4 0.014 18.3 90.0 85.2 
0.025 11.2 84.1 0.015 18.3 91.25 85.1 
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The parameters such equation were determined the method Least 
Squares Mr. Wallace, Division Biology, National Research Council 
Canada. was found that 


where the temperature difference degrees Fahrenheit and the 
turbine output megawatts. The plot this equation shown Fig. 


Temperature Difference,F. 


Turbine Output, Megawatts 


Fic. Turbine efficiency curve. 


The maximum efficiency obtained when 
15,290 kw. 


max. 
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that is, gate opening 65%, Fig. Substituting this value 
Equation (1), and hence the maximum turbine efficiency 
92%. 

The flow water through turbine can determined the temperature 
difference, the height water, and turbine output are known. 


Let horsepower obtainable for height water, ft., 


turbine output, horsepower, 

flow water, c.f.s., 


The temperature difference measured will the temperature difference 
minus the loss temperature difference caused the turbine output. 778 
But 
(2) 


the turbine output measured horsepower, Equation (1) becomes 


From Equation (2) 


Substituting this value for Equation (3) and putting 84.7, obtain, 


The plot this equation shown Fig. and the calculated values 
Table II. 


The most efficient flow 1900 c.f.s. 
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26000 


16000 —<— Most Efficient Flow 


14000 


12000 


10000 


Turbine Output, Horsepower 


8000 


6000 


4000 


2000 


500 1000 1500 2000 2500 3000 3500 
Flow, ¢.f.s. 


Flow water through the turbine cubic feet per second various loads the 
turbine. 


Three curves are shown Fig. that are interest. Curve the plot 
the gate opening against the temperature difference degrees Fahrenheit, 
Curve the gate opening against the turbine output horsepower, and 
Curve the gate opening against calculated flow cubic feet per second. 
The variation efficiency the unit determined this method shown 
Curve The efficiency increases rapidly gate opening about 50%, 
then not rapidly maximum efficiency gate opening 65%, and 


' 
24000 
22000 
20000 
18000 
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Gate opening, 
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Turbine output, hp. 
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TABLE 


DIFFERENT LOADS ON THE TURBINE 


opening about 70%. 


Unit discharge, c.f.s. 


Temperature 


096 
085 
059 
058 
048 
046 
034 
032 
026 
024 


difference, 


TABLE III 


TEMPERATURE DIFFERENCE AND TURBINE OUTPUT AVERAGED FOR EACH SETTING THE GATE 
OPENING, FULL GATE 100%; UNIT DISCHARGE OBTAINED FROM FIG. 


10000 
11000 
12000 
13000 
14000 
16000 
18000 
20000 
22000 
24000 
26000 


Turbine output, hp. 


Turbine output 


Megawatts 


wn 


co 

w 
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CALCULATED UNIT DISCHARGE CUBIC FEET PER SECOND FOR 


Unit discharge, c.f.s. 


1516 
1544 
1604 
1701 
1768 
1912 
2100 
2263 
2483 
2747 
3079 


then drops off the plot obtained from the figures supplied 
the powerhouse operator and indicates maximum efficiency gate 
The values represented dots Curve were 
obtained from Equation (4), the crosses representing the values obtained from 
Equation (2) when the temperature and turbine output are averaged for each 
gate opening, Table 


Unit discharge, 


Horsepower 
1475 780 
2211 950 
4288 1175 
5494 1250 
6700 1325 
7906 1400 
8978 1460 

10318 1550 
11122 
12328 1675 
13132 1725 
14418 1800 
15142 1850 
16147 1900 
16616 1950 
17554 2025 
18760 2125 
20100 2275 
21172 2375 
21172 2375 
22512 2550 
23048 2625 
23584 2700 
23718 2725 
24016 2750 
24321 2800 
24455 2830 
24522 2850 


212 

483 

1000 652 

1400 769 

2000 894 

3000 1104 

4000 1129 

5000 1207 

7000 1339 

8000 1399 

9000 1457 

019 

025 

023 

023 

015 

015 

014 

008 
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005 

009 

009 
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Temperature Curve 
5000 10000 15000 20000 25000 

Turbine Output, Horsepower 
500 1000 2000 2500 3000 
Flow, 


475 


PerCent Gate Opening 


Curve B 


Curve C 


Fic. gate opening the temperature difference between the headrace and the 
A), the turbine output (Curve B), the flow water through the turbine 
urve C). 


90 4 
87-5 
85 
82-5 
80 
77-5 
725 
7O; 4 / 
/ 
| / 
60 4 x 
575 } 
55 ° 
50 a | 
4 
425 
40 a . 
35 / a 
325 O 
30 } 
225] 4 
| 


522 CANADIAN JOURNAL RESEARCH. VOL. 26, SEC. 


was noted that vibration the power unit would begin gate opening 
about 40%, increase maximum about 55%, and then decrease until 
could not felt about 70%. This vibration would cause fluctuations 
the temperature the water the tailrace and shown the scatter 
the values Fig. between the points marked and and also 
the crosses Curve Fig. 
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CANADIAN AERIAL MAGNETIC SURVEYS (M.A.D.)! 


BAILEY 


Abstract 


Tests concluded with airborne magnetic detector measuring the earth’s total 
field indicate high degree accuracy and speed for reconnaissance work over 
large areas. Cheapness operation and freedom from small local anomalies 
make this valuable tool for geophysical prospecting. 


Introduction 


For many years the measurement magnetic anomalies the terrestrial 
field has been valuable method geological studies and mineral location. 
able map field variations rapidly means airborne equipment 
important advance the application geophysics prospecting. This 
now possible result intensive effort during the war develop 
airborne magnetic detector for submarine detection. 


The purpose this report describe the use this instrument. for 
studying magnetic contours over areas known contain magnetic minerals, 
and the method which these data may presented that they may 
useful those familiar with the interpretation magnetic maps relation 
the geology the area. 


Many types instruments depending magnetic field variations were 
tried for detecting purposes during the war. Co-operation between the allied 
research laboratories was always possible and parallel development was carried 
out meet specific requirements the different armed services. Several 
types magnetometers were made the National Research Laboratories, 
Ottawa, and used problems relating harbor protection, automatic boat 
steering, and location lost equipment. 

Experience this field indicated that airborne type would have immense 
value geophysical prospecting; consequently, immediately after the war, 
work was started investigate the best way use for that purpose. The 
most advanced airborne type was that developed the U.S. Naval Ordnance 
Laboratory. Though early experiments were carried out here with magnetom- 
eters designed the National Research Laboratories, was soon apparent 
that the applications these instruments geophysical prospecting would 
proceed most expeditiously obtaining instruments already well developed. 
Accordingly, visits were made Washington that resulted the assistance 
and co-operation available there putting this wartime development into 
practical use. Magnetic detectors were made available the National 
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Research Laboratories Canada for government supported research that 
would further basic knowledge geophysical and terrestrial magnetic 
problems. 


The instruments used the experiment described here were those developed 
the Bell Telephone Laboratories under contract with the U.S. Naval 
Ordnance Laboratory and manufactured the Western Electric Company 
for the Navy Bureau Aeronautics. These instruments were obtained 
through the kind co-operation U.S. Naval Ordnance Laboratory with 
whom friendly liaison had been maintained various other magnetic problems 
during the recent war. our desire here express thanks the U.S. 
Navy, particularly the Bureau Aeronautics and the Ordnance Laboratory, 
for making the Magnetic Airborne Detector (M.A.D.) equipment available 
us. Thanks are also due the United States Geological Survey for their 
kindness sending the M.A.D. equipped plane Canada demonstrate 
this new method survey and for their help introducing survey 
methods and map-making technique. Personal thanks are due Dr. 
Rumbaugh the U.S. Naval Ordnance Laboratories, Commander 
Joyce the U.S. Naval Bureau Aeronautics, Mr. Homer Jensen, formerly 
the U.S. Naval Ordnance Laboratories, and Mr. James Balsley, Jr., 
and his associates the United States Geological 


Parallel with the experiments described here the development M.A.D. 
prospecting instrument was progressing much more rapidly the United 
States. the time writing, hundreds thousands square miles 
country have been covered, and more types instruments are becoming 
available, that their use prospecting already commercial proposition. 


The 


The purpose this instrument detect anomalies the earth’s magnetic 
field associated with the presence oil minerals. The sensitive part 
the instrument consists three magnetometer elements (of the second 
harmonic fluxgate type) mutually perpendicular one another Two 
these elements, operating normal the earth’s field, control servo systems 
which maintain the third element the earth’s total field. The whole 
detecting mechanism towed bomb-shaped housing sufficiently far from 
the aircraft escape the effects the magnetic field the aircraft itself. 
The variations the output the element maintained the direction the 
total field are recorded simple recording milliammeter. 


follow the description refer Fig. Two gimbals whose axes 
rotation are normal each other are driven independently two servo 
motors. are maintained within three-quarters degree the 
direction null field and can corrected with maximum angular velocity 
20° per sec. The inner gimbal carries Lucite cylinder containing the 
three separate magnetometer elements right angles each other. The 
middle element used detector. outside element perpendicular 
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the axis mounting called the transverse element, while the other, mounted 
along the axis rotation, referred the axial coil. When the axial and 
transverse elements are zero field the detector element total field. 
Any movement that destroys this orientation causes the axial and transverse 
elements produce signal proportional the sine the angular displace- 
ment, while the detector element (in total field) produces change signal 
less magnitude since varies the cosine the displacement from total 
field. mechanism, therefore, operating keep the detector element 
total field, activated signal that large comparison the detector 
increment—the one that used for measurement. 


The detector element consists strip cylinder permalloy that 
excited sinusoidal alternatirig voltage carried coil wire around the 
strip permalloy. steady magnetic field the current wave the coil 
resulting from the sinusoidal exciting voltage becomes unsymmetrical and 
prominent second harmonic appears. This harmonic roughly linear 
function the steady field and hence can used for measurement. The 
second harmonic (2000 cycles) separated from the exciting voltage (1000 
cycles) feeding through band pass filter high quality, then the 
amplifier, which serves bring usable level. 

1000 cycle oscillator (V7 and V8, Figs. and set cycles provides 
exciting current through 1000 cycle filters each the three fluxgate 
elements. also feeds doubling circuit (V2) produce 2000 cycles for the 
servo motor’s power. second oscillator (V1) provides output 2130 
cycles cycles), which, when mixed with the 2000 cycles from the 
doubling circuit, provides 130 cycles for one phase each the two phase 
servo motors. The output from the detector element applied amplifier 
(V1 and V4, Figs. and through 2000 cycle filter. Since the sensitivity 
the element per gamma the input the amplifier there would 
signal 2.1 the grid the first tube the total field was 60,000 
gammas. eliminate the difficulty handling large output signal, 
large portion the average earth’s field balanced out direct current 
through the coil the magnetometer element. the present equipment all 
but 1500 gammas neutralized. The remaining average signal 52.5 
mv. amplified conventional a-c. amplifier (V1 and V4) and rectified 
diode (V7), thus producing d-c. potential (across R19). 

the recording sensitivity may adjusted high 100 
gammas for full scale, the amplified signal corresponding 1500 gammas plus 
minus the variations observed still too high feed directly the 
recorder. The signal therefore further controlled imposing d-c. 
neutralizing voltage (across R18) that the output zero for arbitrary 
fixed steady field and produces plus and minus variations from this balance. 
These variations are fed d-c. amplifier (V9 and V10) after going through 
sensitivity control circuit the control box. The output (V10) recorded 
directly. 
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The output from the axial fluxgate passes through 2000 cycle filter 
second amplifier (V2 and V8, Figs. and 7). then mixed with 2130 
cycles (V10 and V5, Figs. and producing 130 cycles, which amplified 
and fed the second winding the servo motor, which acts bring this 
element into zero field. The output the transverse fluxgate passes 
third amplifier and similarly the second winding the other servo. The 
2000 cycle amplified outputs the axial and transverse elements are also 
fed separate rectifiers (V5 and V6, Figs. and 7), yielding d-c. component 
proportional the square the second harmonic signal voltage. These two 
d-c. components are then mixed and fed back through the control box and 
added the detector output appropriate way compensate for small 
and rapid misorientations. 


The control box has been mentioned above connection with the adjusting 
the sensitivity and balance. Its functions are follows. attenuates 
the input the d-c. amplifier that the recorder scales will read 100, 200, 
500, 1500, 5000 gammas full scale. effects the compensation required 
for the small and rapid misorientations referred the previous paragraph. 
provides for the passing d-c. increments through the detector coil 
50-, 500-, and 5000-gamma steps. This accomplished three parallel 
shunts across battery large cell construction. Voltage 
stabilized rapid discharge bring the battery its plateau level. The 
control box load 0.78 mill then permanently applied. The battery 
remains useful for several months. The 50-gamma steps are obtained the 
selection proper fixed points across voltage divider. The 500- and 5000- 
gamma steps are obtained switching d-c. increments through the detector 
coil such manner that the impedance looking into the battery constant. 


The power pack operated from d-c. mains and consists two 
dynamotors parallel, feeding electronically stabilized output circuit. 


Method Adjustment 

Results with this instrument are obtained the following manner. The 
instrument first put operation magnetically quiet location with 
neutralizing d.c. the coil the detecting element. After the servo motors 
place the detector coil total field, direct current gradually increased, 
the detector until the total earth’s field exceeds the neutralizing field 
1500 gammas, i.e., the instrument undernulled 1500 gammas. (This 
1500 gammas arbitrary figure and was chosen since suitable value 
for the type anomalies encountered.) The recorder now reads zero for 
reasons shown above. 


For calibration the instrument set operation and set Helmholtz 
coils placed coaxially with the detector element that known fields may 
added subtracted from the earth’s total field. 

Values indicated the recorder are therefore departures from average 
total field, and the term used this report, refers the graphs 
these variations. 


| 
| 
~ 


BAILEY: CANADIAN AERIAL MAGNETIC SURVEYS (M.A.D.) §27 


The sensitivity the instrument may made correspond 100 gammas 
for full scale the Esterline Angus recorder, various steps 5000 
gammas full scale. 


The absolute value the field, course, not known recorded unless 
specific procedure adopted whereby the instrument calibrated base 
station whose absolute value known. For magnetic survey work this not 
necessary, only anomalies are measured. The compensating voltage may 
maintained quite uniform and the gain the amplifiers reasonably 
stable that the total field probably indicated with consistency 
from day day. Variation over short time can read few gammas 
and drifts are very small when temperatures are steady. 


Note:—Just before the work the airport was terminated, the apparatus 
was operated continuously for few days effort observe diurnal 
changes. The detector was placed weatherproof box 200 ft. from the 
hangar while other equipment remained inside. Temperature has consider- 
able effect the indicated field reading. However, since compensation 
made for such changes the Survey Procedure further time has been given 
the matter. 


Survey Procedure 


The area selected for the experimental survey lay some miles west 
Ottawa. The survey was conducted flying half-mile grids each 
three elevations, 5000, 1000, and 500 ft. The method survey was 
follows. grid lines (north-south) were drawn one-mile-to-the- 
inch topographical maps. (Suitable aerial mosaics were not available.) 


The crew the Anson plane carrying the apparatus consisted the 
navigator, magnetometer operator, and usually fourth member who oper- 
ated the winch, etc. The pilot was responsible for following the grid lines and 
maintaining constant altitude and speed. The navigator traced the plane’s 
course the topographical map and aided the pilot his duty following 
the grid lines indicating land marks. 


While flying the north-south grid, note was made the magnetically 
flat and the magnetically rough areas. This enabled the selection area 
suitable for flying east-west base line bring all north-south lines 
common denominator necessitated diurnal and other magnetic changes 
and instrument drift. 


Sonne continuous strip camera (Fig. was operation during all flights. 
When passing over easily identifiable points, the navigator operated foot 
pedal that linked the magnetic reading with that particular geographical 
position. This was accomplished means photographic and electro- 
magnetic system. 

The recorder used was Esterline Angus. While was satisfactory 
other respects, found the chart rather narrow for the purpose. During 
flight, therefore necessary observe the meter and maintain proper 
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placement the pen applying and recording known increments the 
neutralizing circuit. This results discontinuous profile the recorder 
tape The control box and linkages from the recorder (Fig. operate 
that increments the neutralizing field and changes sensitivity are printed 
the tape. 

reference the developed film, the flight lines are redrawn coincide 
with the actual path the plane. The geographical check numbers are also 
plotted these flight lines from the film. The Esterline Angus tape with 
the magnetic record then placed the graph rectifier This 
machine with system levers and two movable tables convert the 
curvilinear record the recorder tape one with rectilinear co-ordinates, 
(b) expand compress the horizontal scale fit the desired map scale, 
(c) increase decrease the vertical scale, and (d) fit the offset segments 
the magnetic profile the Esterline Angus tape into continuous profile. 
The base line profile values are then established where they geo- 
graphically coincide with the north-south flight line profile values, and the 
comparison noted. there has been terrestrial magnetic change, 
there has been instrumental drift, and the geographical positioning 
exact—then these values will equal. they are unequal necessary 
adjust the reference line for each profile compensate for this 
difference. Flight lines flown over period hours, days, months can 
therefore tied together base line that can flown few minutes, 
thus eliminating the possible errors due drift, temperature, etc. 


Suitable magnetic intensity values are next read from the north-south 
profiles (Fig. 10) and transferred the flight lines (Fig. 11) obtained from 
the continuous strip film. The magnetic highs and lows are noted and 
magnetic contours drawn. 

Three factors must taken into account before magnetic values can 
assigned geographical position. They are distance bird behind camera, 
lag due recording, and perpendicularity the camera.* Probably the 
simplest method for compensating for all three factors fly the base line 
profiles should coincide when tied geographical positions. they not 
coincide, correction must applied all profiles, with due respect 
direction flight, and will equal half the distance between the baseline 
peaks. 


Conclusions 


The final contours are shown Figs. 12, 13, and three levels were 
flown half-mile grids. the working the results this grid seemed 
quite suitable for the 5000 ft. level but was quite inadequate for the 500 ft. 
level. 

The area was flown with the use one-inch-to-the-mile topographical maps, 
which proved wholly unsuitable for wooded areas. sparsely settled regions 


The Sonne camera gyro controlle roll but has compensation for pitch. 
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good aerial mosaics are essential precision results are desired. The com- 
pleted maps can more accurate than the positioning methods. 
would desirable the present strip camera were stabilized for pitch well 
for roll, especially rugged areas where the plane must climb rapidly 
occasions. 

recording radio altimeter AN/APN1 was used, but was not considered 
useful attempt any correlation from the records when the 5000, 1000, and 
500 ft. readings were already plotted. something better than just general 
reconnaissance required magnetic mapping then altimeter much 
greater accuracy would improvement. fact sufficient precision 
were acquired would possible construct topographical maps well. 

While the ‘‘winch, cable, and bird’’ has proved very satisfactory adds 
the need for greater precautionary measures. this respect strongly 
recommended that (a) the amount tow cable released from the winch 
varied—say foot per month flying—so that the flexing occasioned 
the point suspension will not cause fatigue breakage (this flexing 
the bird itself taken universal joint); the whole apparatus 
checked periodically; (c) permanent, positive, semiautomatic device 
installed for cable release case emergency; and (d) the whole system 
subjected periodic 500 static load tests.. 

the time writing, the apparatus has been transferred inboard 
installation Canso aircraft. This eliminates the hazards associated with 
the towed bird system, which used when only small planes are available. 

This new method prospecting places powerful tool the hands the 
mining industry. The slow, laborious, and expensive method point 
point measurement can now replaced cheaper one yielding con- 
tinuous magnetic profile the rate 120 m.p.h. Small surface anomalies 
that might otherwise mask distort ground measurements are lessened, 
completely wiped out, leaving only the general picture from which the over-all 
geological interpretation can made. Drift and diurnal variations are 
minimized the speed operation and the method magnetic baseline 
corrections. There can doubt that the era the dip needle for purposes 
reconnaissance has disappeared, and will future reserved for necessary 
detailed surface work. The way open for complete magnetic world coverage, 
which time will aid the better utilization natural resources. 

modified instrument this type will have extreme value solving future 
navigational problems. Other instruments will designed with sensitivity 
high 0.1 gamma. Further, the principle will applied the con- 
struction instruments for medical purposes, law enforcement, the location 
lost objects, and laboratory and commercial use. The design absolute 
instrument for the recording and measurement the earth’s terrestrial field 
progress this laboratory. 
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Fic. Magnetometer head, showing orienting gimbals and slip ring contacts. The three 
separate magnetometer elements are housed the central Lucite cylinder. 
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CORRECTED FLIGHT LINE BASE 
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THE FUNDAMENTAL ACTION INTENSIFYING SCREENS 
GAMMA 


Abstract 


The nature the intensifying action lead screens and Patterson screens 
has been investigated for gamma radiography. For the lead screens used 
was found that nearly 52% the blackening the film due electrons 
ejected from the front screen and nearly 48% due electrons ejected from the 
back screen. The blackening the film due the absorption the film itself 
the radiation small. Similar results were obtained with Patterson screens 
except that about 10% the effect the screens due low energy radiation 
(light) which was stopped mgm. per cm.? aluminum. 


Introduction 


standard practice gamma radiography place the film between 
sheets lead when making exposure. The sheet lead front the 
film usually has thickness between 0.003 in. and 0.005 in., while the back 
screen may vary from 0.020 in. 0.040 in. thickness. With these screens 
place, the required for photograph given casting very 
much less than required when the lead sheets are not used. For this 
reason the lead sheets are referred intensifying screens. 


The cause the intensification obtained with these screens was investigated 
totally removing the effect each screen turn, and finally the effects 
both them, order find that the film alone. The effect the 
front screen was eliminated separating the screen from the film and then 
applying magnetic field the region between the screen and the film strong 
enough prevent any electrons ejected from the screen from reaching the 
film. For the back screen effect, evacuated space was provided, instead 
back screen, sufficiently long assure negligible effect due backward- 
moving electrons and radiation from the walls the 
combination both methods, the direct effect the radiation the film 
itself was determined. 


Experimental Arrangement 


The source gamma rays used was capsule containing 203.5 mgm. 
radium, which was placed brass container having walls in. thick. 
keep scattering minimum, lead diaphragm, giving total conical angle 
was used, and all interior surfaces the apparatus were lined with light 
cardboard stop scattered electrons without appreciable production 
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X-rays. The source-to-film distance was kept constant in., and all 
exposures were made through in. brass rather than steel since the 
apparatus had nonmagnetic material. The radiographic properties 
brass are close enough those steel make the results applicable 
either. The front lead screen was 0.003 in. thick, and the back screen 0.020 
in. For the Patterson screens the industrial combination No. 245, with the 
thin screen front, was used. Kodak No-Screen film was used throughout 
and developed Kodalk X-Ray Developer, 20° C., for five minutes, with 
agitation the end every minute. The film densities were read 
Marshall Transmission Densitometer. 


provide constant check the work, control exposure was made 
and developed with each set exposures. For this, cardboard film holder, 
having set lead screens, Patterson screens, the same those used 
the other parts the experiment, was placed under in. brass plate and 
given constant exposure with source film distance 9in. constant 
densities obtained for the control experiments indicated consistency 
exposure, development, and film characteristics. well this, unexposed 
strip film was developed with each set exposures give the background 
density. The latter remained constant 0.32. 


Fig. shows the exposure factor versus film density for the 
lead screens (Curve A), and the Patterson screens (Curve with Kodak 
No-Screen Film. These exposures were made with the film firm contact 


with both the front and the back screen, the front screen being fastened 


3.0 


Density 
on 


1.0 


0.5 


| 


Exposure factor 
Fic. versus exposure factor for film. Curve was taken using 
lead front screen 0.003 in. thick and lead back screen 0.020 thick. Curve was taken 


using Patterson Industrial Combination No. 245. The film was exposed through in. 
brass. (B—brass, F.S.—front screen, F—film, B.S.—back screen. 
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brass plate in. thickness. the following experiments, exposure 
factor was used with the lead screens, and with the Patterson screens. 
These exposure factors give density 2.21 each case (see Fig. 1.) 


Part 


Lead Screens 


(a) The Effect the Front Screen 


Consider magnetic field applied parallel the front screen and perpen- 
dicular the incident beam gamma rays, Fig. All ejected electrons 
will move circular paths the magnetic field, the curvatures the paths 


Screen 


Back Screen 


Fic. The schematic arrangement the apparatus used for determining the effect the 
front 


being dependent the energy the electrons and the strength the 
magnetic field, and their excursions from the front screen the direction 
ejection from the screen into the field. Those that are ejected tangentially 
the front screen towards the right (path Fig. will the hardest 
prevent from reaching the film that follow path having radius 
curvature greater than and that are ejected tangentially (path will 
strike the film. the distance between the film and the screen and 
had the value 2.35 cm. the apparatus used. follows then that electrons 
with energies corresponding values less than (where the flux density) 
would all prevented from reaching the film regardless their directions 
emission. Some electrons (see Fig. will ejected backward 
direction from the film itself, and these will curved back into the film 
the magnetic field and produce some darkening. 


Fig. shows the experimental arrangement used determine the effect 
the front screen. and are the poles electromagnet capable 
giving field 5400 gauss across gap 3cm. The film holder consisted 
air-tight metal box, which could evacuated cut down the scattering 
from the air the box. The front screen was fastened the in. brass top 
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the box, while the film and back screen were supported series spacers. 
These spacers made possible vary the separation between the front 
screen and the film from 1.50 cm. 6.50 cm. series five steps. 


Gamma Rays 


Screen 


Fic. Diagram the apparatus used determine the effect the front screen. 


Some the exposures made determine the effect the front lead screen 
are listed Table Only the average values the densities obtained are 


TABLE 


SUMMARY OF RESULTS OF THE DETERMINATION OF THE EFFECT OF THE LEAD FRONT SCREEN 


Separation between 
front screen and the Magnetic field, gauss Bp, gauss-cm. Average density 
film, cm. 


2.35 1.21 

2.35 1500 1760 1.05 

2.35 2790 3280 0.93 

2.35 4240 4980 0.89 

5385 6325 0.89 


given. The individual measurements agreed within about 5%. These 
results appear graphically Fig. 

seen from Table that for zero magnetic field, change the distance 
between the front screen and the film from cm. 2.35 cm. causes decrease 
density from 2.21 1.21. considerable decrease expected 
the distance increased. Also will seen from Fig. that the density 
constant for equal greater than 4980 gauss-cm. 

correct the density maximum magnetic field zero centimeters 
separation, series exposures was made with separations between the front 
screen and the film 1.50, 2.35, 3.70, 5.00, and 6.50 cm. using the apparatus 
shown Fig. all these exposures, the maximum field 5385 gauss was 


Brass 
Screen 
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Density 


2000 4000 6000 8000 
Bp, gauss-cm. 


Fic. Graph showing the variation film density with magnetic field when the electrons 
ejected from the front lead screen are curved back the application the magnetic field shown. 


applied the region between the front screen and the film. The result 
this series exposures shown Fig. which density plotted against 


Density 


Separation, cm. 


Fic. Graph showing the variation the film density with the separation the 
front lead screen and the film when all the electrons ejected from the front lead screen are 
prevented from reaching the film the magnetic field. 
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the separation between the front screen and the film. For all the points 
plotted the magnetic field 5385 gauss.was sufficient ensure the return 
all the electrons emitted the front screen the forward direction. 
That is, these points correspond points beyond the bend the curve shown 
Fig. From Fig. seen extrapolation that for zero separation 
the density would 1.08. This means that removing the effect the front 
screen decreases the density from 2.21 1.08, reduction 51.2%. 


From that portion the curve Fig. where the density constant 
possible determine the maximum energy the electrons emitted from 
the front screen owing the action the incident gamma rays. The value 
where the density becomes constant 4980 gauss-cm., which corresponds 
electron energy 1.07 Mev. for electron ejected tangentially 
2.5 Mev. for electron projected normally. The maximum energy the 
incident gamma ray from radium 2.2 Mev., which can eject from lead 
photoelectron energy 2.1 Mev. The probability the photoelectron’s 
being ejected the forward direction considerably greater than angle 
one considers that the 2.2 Mev. component the gamma 
rays from radium very small, see why the bend the curve Fig. 
appears the region indicated. 


(b) The Effect the Back Screen 


was impossible use magnetic method determine the back screen 
effect owing the fact that there are great many more electrons ejected 
from the film the forward direction than the backward direction. These 
would curved back into the film the magnetic field, and cause very 
great increase the density the film. trial experiment the density 
obtained with field was 1.42 and when the maximum field 5385 gauss was 
applied the density rose 1.62. For this reason, was decided remove 
the effect the back screen the use distance only. The apparatus 
represented Fig. was used. included bell jar, in. diameter and 
in. high, sealed onto in. brass plate. The front screen was fastened 


the brass plate and the film held firm contact with it. The bell jar was 
evacuated for all exposures, with all other conditions the same for the front 
screen determination. 


Exposures were made also with the back screen contact with the film 
serve check, and densities 2.21 were obtained. This value 
agreement with Table and Fig. for exposure factor 55. Then the 
back screen was removed entirely, leaving only the film and the front screen, 
and several additional exposures were made. These gave densities that 
averaged 1.13. This means that the removal the back screen reduced 
the density from 2.21 1.13, reduction 49.0%. 


O00 
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Stopper Vacuum Pump 


-Bell Jar 


Gamma Rays 


the apparatus used determine the effect the back screen. 


(c) The Effect Film Alone 


determine the effect the film alone, combination the two methods 
was employed. Fig. shows the apparatus employed; made use 


Rays 


Iron For 
Magnetic 
Shielding 


Fic. the apparatus used for the determination the effect the film alone. 


air-tight box in. high having brass top in. thick which the front 
screen was attached. The film was supported in. below the front screen, 
and magnetic field 5385 gauss was applied the region indicated. 
heavy iron bar was placed around the box the region the film act 
magnetic shield and thus provide region zero field the position 


Front 
Screen- 
Brass 
Front Screen 
Brass 
Pump 
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the film. necessary prevent electrons ejected from the film from 
curving back into and thus causing increase density. The bottom 
the box was in. below the film, and back screen was used. The box was 
evacuated for all exposures. 


Several exposures were made with the same conditions the other parts 
the experiment described above. These gave average density 0.08 
when, all cases, the background 0.32 was subtracted. When the 
effects both the front and back screen are removed, the density falls from 
2.21 this indicates that the percentage effect the film alone 3.6%. 


Part 


Patterson X-Ray Intensifying Screens, 
Industrial Combination No. 245 


(a) Front Screen 


The effect the front screen was determined the same manner for the 
lead front screen but with the addition one further step. This consisted 
taking additional series exposures, varying the magnetic field before, 
but having the front screen covered with thin sheet aluminum foil (3.05 
mgm. per This thickness foil would not absorb any appreciable 
number the electrons from the front screen yet would absorb the light 
emitted these screens when they are irradiated with gamma rays. Fig. 
shows the results both these series exposures, Curve with the front 
screen alone, and Curve with the front screen covered with aluminum foil. 


Density 


Bp, gauss-cm. 


showing the variation film density with magnetic field when the electrons 
ejected from the front Patterson screen are curved back the application the magnetic field 
shown. For Graph the Patterson front screen was covered aluminum foil 3.05 mgm. 
per thick. For Curve aluminum was used. 


| 
1.1 
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Fig. shows the result varying the spacing between the front screen and 
film, and analogous Fig. for the lead screens. From this graph 
seen that density 1.05 would obtained cm. separation all the 


Density 


Separation, cm. 


Fic. Graph showing the variation the film density with the separation between the 
front Patterson screen and the film when all the electrons ejected from the front screen are 
prevented from reaching the film the magnetic field, and radiation from the front screen 
filtered 3.05 mgm. per aluminum. 


electrons from the front screen were prevented from reaching the film and 
radiation from the front screen that can absorbed 3.05 mgm. per 
aluminum were stopped. This corresponds reduction density 
1.15 and gives Patterson front screen effect 52.3%., From the displace- 
ment the two curves Fig. seen that the light energy that can 
stopped 3.05 mgm. per cm.? aluminum makes 9.8% the total 
effect the Patterson front screen. 


(b) Patterson Back Screen 


The effect the Patterson back screen was determined exactly the same 
way for the lead back screen. Exposures with both front and back screens 
present gave density 2.20, and, when the back screen was removed and 

chamber employed, the density fell 1.15. This means reduction 
density 1.05, and gives Patterson back screen effect 47.8%. 
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(c) The Film Alone 


The effect the film alone was determined also for the Patterson screens 
order obtain check the work. The determination was made 
exactly the same manner when using the lead screens except this case 
the front screen was covered with aluminum foil. The density obtained 
using exposure factor was 0.075; this gave effect 3.4% for 
the film alone. 


determine whether bombardment the lead screen gamma radiation 
produces any appreciable amount radiation the visible the ultra- 
violet regions, two further exposures were made with the lead screens. The 
first one was made with separation 2.35 cm. between the front screen 
and the film with field 5385 gauss and the second similarly but with 
3.05 mgm. per aluminum contact with the lead front screen. These 
two exposures gave the same density; this indicated that any radiation 
was produced the lead front screen, was penetrating enough pass 
through 3.05 mgm. per aluminum. 


Discussion Errors 


The gross percentage error for the experimental work was about 4%, 
which made several factors. source-to-film distance could 
1/16 in., and the exposure times (21.9 min. for the lead 
screens, min. for the Patterson screens) could reproduced sec. 

The processing the exposed film was under very close control, shown 
the fact that the density the control exposure and the unexposed back- 
ground strip never varied more than 0.02. 


Scattered radiation would increase the film density all cases over that 
due any direct radiation and the electrons themselves. However, with 
the collimator used all parts the experiment, and with all inner sur- 
faces the chamber lined with cardboard and the chamber evacuated, the 
scattered radiation was certainly reduced low value. 

the front screen determination the film was situated strong magnetic 
field that any electrons that were ejected backward from the film were 
curved back into and caused increase density. This effect certainly 
contributes some error but should small since comparatively few electrons 
are ejected backward. the arrangement used determining the effect 
the film alone, the larger size the apparatus made possible provide 
magnetic shielding for the region the film, and thereby largely prevented 
increase density due any curving-back effect. 

Since the film had removed from the front screen finite distance 
for the front screen and the film alone effects, indeterminate error was 
introduced. the case the front screen effect, while was possible 
extrapolate the curves Figs. and cm. separation, proof could 
obtained, owing the lack sufficiently powerful magnetic field, that 
the curve continued was assumed. 


‘ 
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further error introduced calculating the percentage effects though 
the density versus exposure factor curves were linear, whereas they are 
slightly curved the region used. 


Discussion and Summary Results 


Table gives summary the effects the front and the back 
screens, and the film alone, for both lead screens and Patterson screens. 
The fact that the three effects add more than 100% itself indicates 
the magnitude the joint errors involved. 


TABLE 


SUMMARY OF ALL THE RESULTS FOR LEAD SCREENS AND PATTERSON SCREENS 


screens 


Lead screens 

| 

| 


Exposure factor used 
Maximum density 2.20 
Front screen effect: 

(a) Due electrons, 42.5 

(b) Due light energy, 0.0 9.8 

(c) Total effect, 51.2 52.3 
Back screen effect, 49.0 47.8 
Film alone effect, 3.6 3.4 

| 
Total 103.8 103.5 


check the percentage the Patterson back screens effect due 
light energy could obtained. However, seems safe assumption say 
that the same order magnitude that for the front screen. 


The effects reported here were obtained for only one thickness lead 
screens and one type Patterson screen, and for one thickness 
probable that these factors were changed the same general type results 
would obtained. 

obvious from this table that the use some form intensifying screen 
imperative gamma radiography. the effect the film alone should 
small understandable when realized that the only which 
rays can cause darkening the ejection electrons, and only 
those having paths largely within the emulsion will produce any 
appreciable effect. Since the absorption the film gamma rays small, 
comparatively few electrons are available for darkening. doubtful 
whether this fact has been fully appreciated film manufacturers users. 
The absorption the gamma rays the film could increased materials 
high atomic number were incorporated the emulsion. The fact that 
most the darkening comes from the effects the front and back screen 
indicative the care that should taken the choice these screens. 


me 
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The practical aspects the choice these screens from the point view 
the sensitivity obtainable and the intensification that can secured has 
been discussed some the authors another publication (1). will 
suffice state here that the screens should maintained firm contact 
with the film that the electrons ejected from the screen will affect the film 
their point origin, rather than considerable distance, avoid loss 
sharpness. the second place the screens should clean. small 
section the screen covered with some foreign material, the latter will 
absorb the electrons from the screen and will produce only smaller number 
electrons owing its own absorption the gamma rays. This will result 
light spot the film. 
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THE EXTENDED USE KLEINLOGEL’S 


Abstract 


This paper contains but single idea. Its purpose show how, the use 
Kleinlogel’s Rahmenformeln, the number simultaneous linear equations 
that occur the analysis indeterminate structures can reduced materially 
using primary structure, itself made indeterminate parts, provided 
the bending moment diagrams can drawn for them. These are obtained 
from the Rahmenformeln. order recall the reader’s mind the well 
known moment-area method brief summary that process included. 
practical numerical example the analysis frame indeterminate the 
seventh degree worked out detail. arrangements the indeterminate 
components are chosen the illustrative example. The final bending moment 
diagram shown for simple loading case and the influence line for the bending 
moment one point discussed briefly. 


Introduction 


During the recent war the Eighth Edition Professor well 
known Rahmenformeln was reprinted the United States the authority 
the Alien Property Custodian. was published Frederick Ungar 
Publishing Co. 1943, supplemented with crib English—published 
separately—for the benefit those not familiar with the German language. 


This remarkable book contains diagrams and formulae for simple frames 
various shapes, support conditions, and loadings. One hundred and four- 
teen frames are listed with 1643 diagrams accompanied numerous formulae. 
The diagrams show the moment-curves for each case loading; the formulae 
give the algebraic expressions for the moments, shears, and reactions. The 
index contents odd one consisting, does, series simple 
illustrations showing the shape and support conditions each frame with 
the numbers the pages which the moment-curves and formulae pertaining 
will found. very brief introduction contains all the necessary 
information concerning notation, rules for signs, and other pertinent inform- 
ation. However, one does not need able read this part make use 
the tables, which are quite obvious meaning. Besides, the English 
translation now available obviates this difficulty. 


first sight, one might think that this useful handbook would limited 
its application only those cases contained it. This is, however, 
false conclusion for, although not indicated the introduction, the 
simple frames can used units with which build multiple frame 
structures unlimited extent. such synthesis, the use and value the 
book becomes considerably extended and the purpose this paper 
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present method which the analysis multiple frame structures may 
carried out simple and direct manner based the information contained 
it. 

This method, namely, the moment-area method, consists dissociating 
the multiple frame into units that can found the text and using these 
units the elements primary structure. The elements themselves are 
statically indeterminate, but, since the analysis for them already hand 
the text, difficulty way preliminary analysis present and the 
synthesis carried out with these statically indeterminate units instead 
the statically determinate units usually used for that purpose. The amount 
effort required analyze highly indeterminate structure enormously 
reduced and for that reason the analysis multiple structural frames, which 
might appear formidable, brought within quite reasonable range. 


The Moment-area Method 


The purpose this section recall the mind the reader the moment- 
area method. 


order demonstrate this most useful, perhaps, all methods struc- 
tural analysis, simple algebraic example will given. Consider rigid frame 
shown Fig. with which the reader will undoubtedly quite familiar. 


Fics. and 


This structure indeterminate the third degree, three unknown, 
redundant, reactions addition the three reactions plane statics are 
present, and, therefore, addition the three equations equilibrium 
plane statics, three additional equations terms the three redundants must 
established. These form system independent, simultaneous, linear 
equations that must eventually solved. The moment-area method one 
that sets these equations; does not offer method solution them. 
When they are few number the solution them can easily and quickly 
accomplished. 


The first step the process render the indeterminate structure deter- 
minate the insertion suitable cuts and hinges. the present case, 
this shown Fig. which three hinges have been inserted. The struc- 
ture now statically determinate and stable. This called the 
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From the equations virtual-work, the set simultaneous equations 
involving the redundants M,, and the unknown moments the 
selected points the indeterminate structure, are readily established. These 
may written condensed form 


Aji jito3 


The are the loading terms. easily shown that the a;; are given the 
“work owing self-stressed condition. The most important 


Case 


+ 


Case 


term the complete work equation and the only one usually retained, but not 
necessarily so, analysis this sort 


0 
this equation, the and the are the bending moments each mem- 
ber due self-stressed condition imposed placing successively each 
hinge couple, The the product the modulus elasticity 
the moment inertia the cross section for each member. Since often 
constant and occurs both sides the equation, will assumed 
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disappear from the following equations. The integration taken over the 
length each member and the summation sign indicates that the sum these 
integrations must taken. 

This process requires three bending moment diagrams, two which are 
shown Figs. (1), and (2), respectively. The third like (1) except 
reversed. The areas displayed the diagrams are called moment-areas 
and the method takes its name from this fact. 

The simple integrations for many cases have been worked out and are 
given table Die Graphische Statik der Baukonstruk- 
tionen, Band II, abt. page 56. Many combinations moment-areas are 
included this table but here only one will given which all that 
necessary for the present purpose. 

Fig. indicates the case referred and shows two trapezoidal moment- 


areas that are multiplied and integrated from The factor 


will considered constant for the member and, therefore, omitted the 
integration. Then 

The moment-areas and may said this process. The 
matrix speak, represents the structure mathematical form and 
might called the depends, however, the choice the 
primary structure. Owing the assumptions Hooke’s law and elasticity 
the members the structure, the matrix necessity always symmetric. 

Returning the analysis, have now work out the six constants 

pairing the moment-areas for the frame. 


sample calculation, take 


Starting with the left leg the frame, 


l 


2 
| 
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and proceeding around the frame, one finds 


the same way, all the other coefficients may set down, practically 
inspection. Thus, have: 


customary indicate the ratio the moment inertia member 
its length single letter, called the Here, however, 


shall which concordance with Kleinlogel’s formulae. 


The loading terms are determined also from the moment-area the primary 
structure, due any specified loading it, paired with each the three self- 
stressed moment-areas. If, for example, the loading shown Fig. 


The equations then become, matrix form, 


q 
| 
| | 
P 

: 3 / 
~ 


MORRISON: EXTENDED USE RAHMENFORMELN 557 


This completes the demonstration, for there point solving these 
algebraic equations. The object has been merely show how they are set 
means the moment-area method. course, much more ex- 
pedient work directly with the numerical values. 


Solution the Equations 


The main difficulty that arises any method analytical analysis 
highly indeterminate structure the solution the large set simul- 
taneous equations that arise the process. This becomes very tedious and 
great care required checking the solution progresses. For this 
reason any method which the number equations can reduced offers 
material advantage. 


There are number methods that have been developed which the 
number simultaneous equations can reduced. The present paper 
presents such method. based, however, the assumption that the 
moment-area diagrams for number simple, though statically indeter- 
minate, frames are readily hand mentioned the introduction. Klein- 
logel’s Rahmenformeln provides this information. 


The Use Kleinlogel’s 


now turn the use Kleinlogel’s diagrams. the preceding moment- 
area method, the primary structure was chosen statically determinate solely 
because was easy construct the bending moment diagrams for such 
structure. the Rahmenformeln, however, have access literally 
hundreds cases simple, indeterminate frames. Instead, then, rendering 
the statically indeterminate structure determinate, obtain tractable 
primary structure, may adopt primary structure made frames the 
moment-areas which are already known and therefore also easily handled. 


This method dissociating the structure into primary structure, itself 
composed indeterminate parts, reduces very materially the number 
simultaneous equations solved, will seen the following examples. 


Numerical Examples 


order illustrate the advantage gained this process, proposed 
analyze the relatively simple structure shown Fig. the 
primary structure decided upon. This may done several ways. 
this case, the choice shown Fig. The original structure indeter- 
minate the seventh degree. inserting two hinges shown, the primary 
structure indeterminate the fourth degree.* selecting the primary 
structure care must exercised select only such arrangement hinges 
cuts will leave the elements the primary structure such that the 
moment-area diagram already known. should noted that the shorten- 


The left-hand hinge horizontal shelf, speak, and removes two constraints. The 
right-hand hinge removes only one constraint. 


7 
q 
4 


558 CANADIAN JOURNAL RESEARCH. VOL. 26, SEC. 


ing the posts assumed neglected usual and, consequence, 
there will vertical displacement the hinges. solution approxi- 
mate this respect. 


H M, 
Fic. 


The values shown Fig. have been obtained first selecting 
one value the stiffness factor unit value and then finding the ratios 
the respective stiffness factors the remaining members. the 


present case the girder the center span was chosen that 


Each pairing the moment-areas multiplied the corresponding k-value 
for the member. 


Since the structure indeterminate the seventh degree would, 
the primary structure were rendered determinate the usual procedure, 
n(n 

terms; itself formidable task. Then seven simultaneous linear equations 
would have solved; this time consuming, say the least. 


have determine independent coefficients, and loading 


using the primary structure composed elements, themselves 
indeterminate, the number equations solved has been reduced three 
and the number coefficients six, thus effecting very large saving 
effort. 

any case, determine the number redundants, let the degree 
indeterminateness the structure and the sum the degrees indeter- 
minateness the components the primary structure, then the number 
redundants, say 


the present case, the left-hand girder indeterminate the first degree*, 
the middle frame the third degree, and the right-hand girder the zero 
degree. Thus, 


Note that the right-hand end roller supported the post. 
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which the number redundants selected and shown Fig. 
and The arrows indicate the positive directions the forces and couples. 

Figs. (1), (2), (3) show the moment-areas for the structure obtained 
from the Rahmenformeln for the self-stressed conditions —1, 
—1,* respectively. pairing the diagrams, already explained, the 


Loading Case. 
Fic. 


numerical coefficients are easily obtained. They are shown elements 
the matrix. 

The equations have been solved for the simple case loading shown 
Fig. The reciprocal matrix for three equations very easily obtained 
and once set down serves for all possible cases loading and temperature 
changes. 


Actually the connecting force, the right-hand end, the left span the center span. 
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0.7927 0.1223 —1.226 
0.1223 0.5115 1.226 
1.226 1.226 352.7 


—44.64 0.7770 
278.1 
.7770 0.3905 


3397.77 


135.3 


849.55 


135.3 6.27 


14.763 


Ok 
135.3 0.109 


The final bending moment diagram for the loading chosen shown Fig. 
The bending moment diagrams Fig. are used determine the bending 


moment diagrams Fig. 


second analysis the same structure which cut has been made 
the middle the center span shown Figs. and 10. The equations 
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Case 


‘ 


+3.48 
Loading Case 


Fic. 10. 


follow and the results shown Fig. check within reason those the first 

analysis. All the numerical work, except for the final reciprocal matrix, was 
s . . 

carried ordinary slide rule. 
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cases symmetrical structures, advantage should taken the 
symmetry that can maintained suitable choice cuts hinges. Also 
the method group loadings can used advantage wherever group 
loadings tend simplify the moment-area diagrams. group loadings 
meant the choice two more unit loads, couples, applied the primary 
structure simultaneously produce self-stressed condition. This method, 
however, does not reduce the number redundants. 

The method the elastic center can applied also making the 
coefficients equal zero, but preliminary investigation shows that, 
except the case symmetrical structures, there little, any, advantage 


Influence Lines 


further detailed study the structure required, influence lines should 
drawn. This can done setting the equations the influence line 
for any chosen item that requires detailed attention. 


illustration, let required draw the influence line for 
the solution, only the loading terms will change; therefore, the equation for 
will be, from the reciprocal matrix, 


which the C’s are determined pairing the moment-area diagram for 
the unit load with each the moment-area diagrams for the three self-stressed 
conditions. 
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The moment-area diagrams for the unit load each span successively are 
shown Fig. will one equation for each span. For Case 
only will present and for Case only For Case however, all 


three values will exist. 


Coses Line 


12. 


Fic. 13. 


order express the position the moving unit load, shall take 
distance from the left-hand end each span, and from the right-hand 
end. For the unit load the left-hand span can most readily obtained 
from Fig. 13. Thus, from pairing the moment-areas shown, find 
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similar way obtain for the center and the right-hand spans respectively 


Inserting these values Equation (1), find for the three spans, 


VOL. 26, SEC. 


and 
—0.1132 xx’ 
which 
7 = a, - = b. 
respectively, 
MM, = 


The plotted influence line 


shown Fig. 14. 
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PRECISION HEAT TRANSFER MEASUREMENTS WITH 
THERMOCOUPLES—INSULATION 


Abstract 


has been developed for calculating the temperature variation 
insulated thermocouple lead wires that not follow isothermal path. 
The difference between the temperature the junction and that the surround- 
ing material that purports measure has been called error.’ 
has been shown that insulation error determined variations the tempera- 
ture the path followed the lead wires only over limited distance from the 
junction, which has been called the eliminate 
insulation error the path the wires need isothermal only for the critical 
distance. simple method has been developed for calculating the critical 
distance and the insulation error. When the path the wires cannot made 
isothermal the conditions for minimum experimental error are shown small 
diameter wires low specific conductivity with minimum insulation. 


The determination fluid film heat transfer coefficients requires the 
measurement the surface temperature the tube other dividing wall 
through which heat being transferred. This usually requires that 
thermocouple embedded the tube wall and the surface temperature 
calculated from the temperature indicated the couple. 

previous paper (1) the effects geometric errors determining the 
location the thermocouple tip with respect the tube surface were analyzed 
and methods minimizing them were described. The other major source 
error such systems may termed since occurs when 
insulated lead wires are carried away from the junction along path that 
not isothermal. then transmitted to, from, the junction along the 
wire itself that the temperature the junction actually different from 
that the material surrounding it. has long been recognized qualitatively 
that the lead wires should brought out along isothermal path. But 
practice not usually possible adhere strictly this condition, parti- 
cularly the case tubes, for which chordal holes are preferably employed. 
hoped that the present quantitative analysis will help investigators, not 
only estimate the magnitude errors that cannot avoided, but more 
particularly choose experimental arrangements that will minimize these 
errors. 


Qualitative Discussion 


convenient discuss the problem first terms single insulated 
wire completely embedded solid metal (e.g., low-melting solder) 
chordal hole tube wall, such that Fig. 

Heat will flowing radially through the tube wall and the temperature 
gradient within the wall will depend the density heat flux and the 
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thermal conductivity the wall. Table shows that the temperature 
gradient the wall may quite steep and that the temperature difference 
across the width fine thermocouple lead wire may amount several 


Metal 


Copper 


Aluminum 


Yellow brass 


Steel 


Thickness 


TIP COUPLE 
EMBEDDED SOLDER 


CENTERING 
BUSHING 


TUBING 


Fic. thermocouple installation. 


TABLE 


TEMPERATURE GRADIENTS THROUGH METAL WALLS 


400° Heat flux, Temp. Temp. difference F.) across 

2605 400 ,000 153 20.4 1.6 

200 ,000 10.2 0.8 

100 ,000 0.4 

2.6 0.2 

1560 400 ,000 256 34.0 

200 ,000 128 17.0 

4.2 0.3 

1010 400 ,000 396 52.6 4.2 

200 ,000 198 26.3 2.1 

100 ,000 1.0 

6.6 0.5 

300 200 ,000 667 88.7 7.0 

100 ,000 334 44.4 


Kent, ‘‘Mechanical Handbook", Vol. 3-28. 
Thickness 0.133 in. 
0.0105 in. 
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degrees. result this temperature difference across the wire, heat 
will flow through the insuiation into the wire, across the wire, and out again 
through the insulation the other side. Because the high thermal 
resistance the electrical insulation, the amount heat that takes this 
path will small relative the flux through the tube wall itself, but may 
not negligible. Furthermore, the wire not isothermal along its length, 
some the heat that enters through the insulation given point will 
flow along the wire before passing out through the insulation the other 
side. The components flow across the wire and along the wire may 
considered separately. The former component not major importance 
the present discussion does not alter the temperature the wire, 
but the component along the wire causes the temperature the wire differ 
from that its immediate surroundings, and the cause insulation error. 


The first conclusions importance will made clear reference 
Figs. and which illustrate wire subjected only transverse component 
heat flux. Fig. shows cross section the tube and wire the plane 


OUTER SURFACE OF TUBE 


DIRECTION OF FLUX 


tSOTHERMS 


ph 


NOTE: THICKNESS INSULATION INNER SURFACE TUBE 
EXAGGERATED 


Fic. cross section tube plane heat flux showing distortion isotherms 
insulated wire—for case longitudinal heat flux. 


the heat flux; this longitudinal section the case the tube; for 
simplicity throughout this paper will assumed right cross section 
the case the wire, although strictly speaking the axis the wire not 
necessarily normal the radii the tube, that cross sections the wire 
the plane the heat flux may actually slightly elliptical shape. 
This figure illustrates roughly how the isotherms, normally parallel, are 
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NOTE - THICKNESS OF 
| ' INSULATION IS EXAGGERATED 


TEMPERATURE > 


TEMPERATURE 
WIRE 

TEMPERATURE 


— — TEMPERATURE VARIATION iN WALL WITHOUT WIRE 


TEMPERATURE VARIATION CAUSED BY WIRE 


Fic. Temperature variation across wire—for case longitudinal heat flux. 


distorted the region the wire. Temperatures along the line Fig. 
line which passes through the center the wire, are plotted against distance 
Fig. From these figures may draw two conclusions: 


(a) The temperature drop largely concentrated the insulation that 
the temperature the wire metal may considered essentially constant 
across any given cross section taken the direction heat flow; this wire 
temperature any given point will Then let the tempera- 
ture that would have obtained the metal wall the same point (correspond- 
ing the center the wire) the wire had not been there. This ‘‘would- 
temperature the average the actual temperatures existing 
practice the two sides the wire (see Fig. 3), and, since thus the 
average temperature which the periphery the wire any cross section 
subjected, will termed the ambient temperature. 


(b) Assuming that the cross section the wire and its insulation 
symmetrical and that the temperature gradient the surrounding metal 
wall uniform, see from Fig. that when there longitudinal 
distribution heat. 


But when the wire not brought out along isothermal path, heat flows 
along the wire well across it, and longer equal Then the 
tip the wire (i.e., the 
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Insulation Temperature Drop 


any given cross section the wire, the temperature the outer surface 
the insulation varies from maximum the side nearest the inside the 
tube minimum the side nearest the outside the tube. Since the 
temperature variation the tube wall linear with distance one neglects 
distortion the isotherms the immediate region the wire, the projection 
the peripheral temperature the diameter will, first approximation, 
also linear with distance (see Fig. and compare with Fig. Fig. 


DIRECTION FLUX 


INSULATIO 


WIRE METAL 


APPROXIMATE PERIPHERAL 


SLOPE = 


Fic. Calculation peripheral temperature. 


distance measured across the diameter the wire from the hotter side, 
the temperature gradient the tube wall,* and the temperature 
any point the periphery the wire. From the geometry Fig. 


follows that 
(2) 
and, 


Since the temperature the wire metal may considered essentially constant 
throughout the cross section, the temperature difference across the insulation 
any point around the periphery 


AT; (4) 


Since the temperature decreases the positive direction negative. 


} 
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This temperature difference not constant around the wire, but varies 
from maximum positive value one side the wire maximum negative 
value the other side. When the positive and negative values are 
symmetrical and equal, the quantities heat entering and leaving the wire 
the cross section are equal and there flow heat along the wire. When 
the positive and negative values are not equal, the effective causing 
net flow heat into the wire the average value AT; taken around the 
periphery the wire; that is, 


Substituting from Equations (4) and (3), integrating, and then substituting 
Equation (2): 
(6) 


That is, the average temperature difference causing net flow heat into, 
out of, the wire any given cross section simply the difference between 
the ambient temperature and the temperature the wire. 


can seen intuitively that Equation (6) not restricted the hypothe- 
tical case which linear function the integration will give the 


same result provided only that symmetrical with respect the 
center the wire. This condition, and hence Equation (6), may assumed 
approximately true any cross section the wire (see dotted curve 
Fig. 4), although this condition rigorously true only the limiting case 
when the fluxes entering and leaving the wire the same cross section are 
equal, when there longitudinal flux and equals illustrated 
Fig. 


Wire Temperature Equation 


order obtain expression defining the temperature the wire metal 
any point, consider increment the wire taken right cross section, 
inches thick, distance from the tip, shown Fig. general, 
the direction increasing temperature the wire metal will the same 
the direction increasing ambient temperature, since the change 
ambient temperature that causes change wire temperature. The ambient 
temperature will some function and for convenience the temperature 
datum will taken zero ambient temperature the tip the wire. 
Consider the flow heat through the cross section into the part the 
wire the left the section the positive direction heat flux, since this 
has been arbitrarily defined Fig. being the direction decreasing wire 
temperature. 


4 
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THERMOCOUPLE 
WIRE 


WIRE TEMP. CURVE 


CURVE 


TEMPERATURE 


Fic. Calculation wire temperature equation for known ambient temperature curve. 


Then according Newton’s law the flow heat across the increment 


given 


where the symbols are defined the Table Nomenclature. 


Since the flow heat out through the surface the increment 


(8) 
the heat lost through the surface and end the wire will 


Since heat not stored any part the wire, may equate (7) and (9), 
and noting that definition obtain 


(10) 


where (10a) 


and dt/ds the temperature gradient the wire 


| 


572 CANADIAN JOURNAL RESEARCH. VOL. 26, SEC. 


Equation (10) defines the curve f(s) giving the temperature the wire 
any point. practice may differ for different regions along the length 
the wire, best keep inside the integral sign. 


Since the heat lost through the end the wire may usually neglected, 
will shown later, Equation (10) may written 


(11) 
0 0 


And for value corresponding maximum minimum the curve 


(12) 


Equation (12) indicates that the wire may considered divided. into 
separate regions the maxima and minima the curve, and that tempera- 
ture effects one region not carry over into the next. also evident 
that the maxima and minima the curve f(s) coincide with the maxima 
and minima the curve F(s) since the variations are due solely 
the variations and hence must always the same sense those 
This generalization helpful plotting the wire temperature curve, and 
makes Equation (12) applicable where defined the maxima and minima 
the ambient temperature curve. 


Equation (11) cannot solved directly determine the wire temperature 
curve, but Equations (11) and (12) provide the criteria which may 
check any wire temperature curve obtained trial and error approximate 
methods. 


Approximate Solution 


approximate solution Equation (10), sufficiently accurate for all 
practical purposes, may obtained converting the approximate 
form 


where the average value for each increment concerned. 


each increment assumed sufficiently small that the variation 
across may considered linear with distance, may write 


| 
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for any increment. Using Equation (14), the second term Equation (13) 
may broken into some its constituent increments, that Equation 
(13) becomes 


Pn As, (14-2 + Ata-1 4- k to At, 


The significance the various subscripts will made clear reference 
Fig. 


be 


wire 
ent 

TEMPERATURE DATUM: 
z 
DISTANCE DATUM: 
” 
z 

4s, Ass 


Fic. illustrating approximate solution wire temperature equation. 


Regrouping, letting As, that the term containing Af, dis- 
appears, and solving for obtain 


| { 
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which 


Equation (16) may used calculate for each increment turn 
except the For example, when and for case which 


5, 


For the first increment, and Equations (13) and (14) may com- 
bined and reduced 


0 : wk 
which 


4k. As; 


(19) 


the correction for heat passing through 
1 


the end area the wire, 


(20) 
When F(s) may considered linear over the first increment, Equation 


(19) reduces 


(21) 


Value 


Table lists thermal conductivities materials commonly used for 
electrical insulation; the conductivity copper wire the same units 
about unity. For purposes illustration, will assume average value 

Assuming No. B&S enamelled wire with diameter 0.0105 in. 
outside the insulation and diameter 0.0100 in. for the bare wire, the 
thickness the insulation may taken 0.00025 in. These values 
readily determined with traveling microscope the case the 


“4 
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TABLE 
THERMAL CONDUCTIVITIES INSULATING MATERIALS 


cal./(sec.) (sq. cm.) (°C./cm.) 


Polyvinyl formal 
chloride-acetate 3.9-4.0 
Rubber 4.3-4.5 
Gutta percha 4.8 
Varnished cloth 
Polyethylene 


Authorities: Miner, ‘‘Insulation Electrical McGraw-Hill, 
Catalogue,” 1946. Physics and Chemistry.” 


original wire insulation. more difficult estimate the thickness 
insulation such Heresite applied over the tip after soldering the thermo- 
couple. 


For such wire, 0.05, Equation (20). 


evident, therefore, that neglecting the term Equation (19) 
(21) will, the present example, introduce error only 5%, which 
much less than the error that would normally encountered estimating 
and for the insulation, that this factor allow for heat entering the 
wire through its end area, may neglected. With equally negligible 
error the corresponding term may neglected the equations 
where occurs. obvious that there may cases, particularly with 
thicker wires, which not negligible, but for simplicity this end effect 
will neglected throughout the rest this paper. Then Equations (1) 
and (21) reduce to: 


Calculation Ambient Temperature 


The method calculating the ambient temperature depends, course, 
upon the geometry the thermocouple arrangement employed. Since the 
case the chordal hole shown Fig. similar many arrangements 
commonly used, the ambient temperature curve has been calculated from the 
tip the couple the thermocouple groove. was assumed that the tip 
the couple 0.500 in. from the surface measured along the axis the 
hole, and that the wire bent down from the axis the hole distance 
0.45 in. from the tip, reaching the bottom the groove 0.55 in. from 
the tip. The bottom the groove itself was assumed relatively isothermal 
over moderate distances, and the flux was taken such produce tempera- 
ture gradient 50° per inch the metal wall the tube. 
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From the geometry Fig. will obvious that 


0.500 
and that 


The calculated values are shown Table III and the resulting ambient 
temperature curve shown Fig. typical ambient temperature 
curve will used illustrate the calculation wire temperature. 


where 


TABLE III 
CALCULATION AMBIENT TEMPERATURE 


Distance from tip in. inches temperature, 


0.68 
1.18 
1.46 
1.53 
1.38 
1.03 
0.47 
—0.28 


2.37 


ooocooooo 
OS 
ASnonone 


0.40 
0.45 


Temperature datum. 
bottom groove, 0.0937 in. deep; 0.6575 0.0937. 


TEMP. VARIATION FROM 


AMBIENT TEMPERATURE 
WIRE 
WIRE 
WIRE 


WIRE ow BOTTOMOF GROOVE 


DISTANCE ALONG WIRE FROM TIP—— INCHES 
Fic. Typical wire temperature curves. 


0.6112 

0.5877 
0.5820 
0.5806 
0.5835 
0.5906 
0.6018 
0.6169 
0.6356 
0.5638 

-1,0 
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Calculation Wire Temperature 


order give better understanding the phenomena involved, the 
wire temperature has been calculated for three typical cases and the calcula- 
tions have been carried several significant figures that the trends will 
clearly shown. The calculation consists stepwise solution, one increment 
time, proceeding along the wire from its tip. 


Case 

Assume that the insulation throughout the length the wire equivalent 
the original insulation, 0.00025 in. thick. 
Then 800 Equation (10a); 


The calculations are tabulated Table and the resulting wire tempera- 
ture curve plotted Fig. For the first increment one uses the simplified 
form Equation (19), which end effect neglected: 


And for the second increment one uses the similarly simplified form 
Equation (17): 


For subsequent increments Equation (16), simplified because con- 
stant, used the form 


The first step the calculation fill the part the table the left 
the double line, which involves integrating the area under the ambient 
temperature curve for each increment. the present work the integration 
was performed the primitive but nonetheless satisfactory method 
counting squares the graph paper under the curve Fig. One then 
proceeds working across the table from left right, completing each line 
before dropping down the next. this systematic procedure followed, 
what can rather confusing calculation becomes quite straightforward. 


the calculations Case may apply the criterion provided 
Equation (12). The first maximum the ambient temperature curve 
occurs 0.191 (approx. 4), and the minimum occurs 0.450 
9). Applying the criterion, for the first maximum: 


0 
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and similarly, for the minimum: 


0 


TABLE 
CALCULATION WIRE TEMPERATURE CURVE 


| 
| 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


OCH KO 


-351 
748 
.393 
.449 
302 
0.949 


Case III identical with Case down the heavy line 0.35, which point the 
increase insulation thickness occurs. 


Case 
22.00 22.00 59.84 
37.82 59.82 91.44 
50.72 110.54 111.44 
56.84 167.38 115.92 
55.02 222.40 104.16 
45.02 267.42 75.92 
26.82 294.24 31.36 
0.66 294.90 28.72 
1.94 277.14 57.84 
50.88 328.02 145.68 

89.92 498.10 184.72 
93.16 591.26 187.92 
94.62 780.12 189.44 
18.17 18.17 46.76 
24.61 42.78 51.68 
20.18 62.96 29.04 
9.14 76.80 46.80 
31.40 108.20 78.80 
42.07 150.27 89.48 
46.81 242.71 94.20 
*Case III 
0.10 1.18 51.88 0.397 22.00 22.00 59.84 
0.15 1.46 105.52 0.395 59.82 91.44 
0.20 1.53 166.24 0.250 50.72 110.54 111.44 
0.25 1.38 225.32 0.056 56.84 167.38 115.92 
0.30 1.03 274.48 —0.147 55.02 222.40 104.16 
0.35 0.47 305.40 —0.353 45.02 
0.982 2.37 325.90 0.958 1.890 8.92 313.96 23.89 
1.298 340.88 0.320 2.210 12.96 326.92 27.93 
1.614 355.86 0.107 14.31 341.23 
1.930 2.37 370.84 0.035 14.76 355.99 
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The figure 1.14 error 0.7% the first maximum, and 5.36 
error the minimum. These errors are primarily due the assump- 
tion which Equation (14) based, namely, that the wire temperature 
curve may considered linear over the length each increment. will 
noted that the error introduced this assumption not large. 
increase would tend increase this error, but when increase 
corresponds increase the amount insulation tends make 
the curve flatter and thereby has compensating effect. 


Case 

Assume that the insulation throughout the length the wire 0.001 in. 

The calculations here are identical with those for Case and are shown 
Table IV. The resulting wire temperature curve shown Fig. indicates 
the effect four-fold increase insulation thickness compared with 
Case 

interesting this point compare the values insulation error 
obtained the approximate equation (22) with those obtained the more 
accurate form (23). For Case the approximate equation gives insulation 
error 1.18/2 0.590° compared with the more accurate value 
0.648° F., error 9%. For Case the corresponding figures are 0.340 
and 0.351, error 3%. obvious, therefore, that the approximate 
form (Equation 22) may only used when relatively short, that is, 
when the insulation quite thin. 


Case III 

Assume the case, illustrated Fig. which the insulation thickness 
0.00025 in. from the tip the wire the back the centering bushing—a 
distance that may taken 0.35 in.; and assume that the air insulation 
the drilled hole back the bushing equivalent insulation thickness 
0.01 in. case similar the installation described the previous 
paper (1) and is, therefore, typical what may encountered practice. 

The first part the calculations identical with that for Case For the 
second section the wire, 20; 0.316; and 6.32. 
These calculations are also tabulated Table IV. the transition from 
one region the other, care must taken use the correct values and 
the equations. (16) must used for calculating and 
must used calculating Otherwise the procedure 
identical with that employed for Cases and 


Critical Distance 


examination Fig. shows that the wire temperature approaches the 
temperature the thermocouple groove more slowly the thicker the insula- 
tion. Similarly, disturbances the ambient temperature (such the 
minimum 0.45) have less effect wire temperature when the insula- 


4 
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tion thicker, and the disturbance wire temperature then distributed 
along greater length the wire. This indicates qualitatively the fact that 
the distance along the wire that disturbance ambient temperature 
transmitted some function the thermal resistance the insulation; 
would also expected function the conductivity the wire. 
Combining Equations (10a) and (16a) gives just such relation, namely, 


(26) 


for w/kD proportional the thermal resistance the insulation per unit 
length wire, and KD? proportional the conductivity the wire per 
unit length. therefore see illustrated what has already been demon- 
strated the mathematics the analysis, namely, that defines the 
maximum distance that disturbance specific point can transmitted 
along the wire. will, therefore, called the 


The significant conclusion therefore, reached that the path the lead wires 
from the junction need isothermal only for the critical distance; and the 
path not quite isothermal, then the insulation error determined the variation 
ambient temperature only over the critical distance. 


For Cases and above, the critical distances are 0.05 in. and 0.10 in., 
which are surprisingly short. obvious from Equation (26) that the 
critical distance, and hence the insulation error, larger, the larger the 
thermal resistance, w/k, the insulation. follows, therefore, that experi- 
menters should strive for minimum thermal insulation the thermo- 
couple wires near the tip without losing the necessary* electrical insulation. 
Under most circumstances unusually thick blob insulation the tip 
thermocouple could seriously affect the accuracy the readings. Other 
factors that will tend reduce insulation error are small wire diameter and 
low specific conductivity wire metal. 


Effect Joining Dissimilar Wires 


The discussion thus far has been terms single wire, whereas thermo- 
couple consists two dissimilar wires, having different values and 
perhaps Applying Equation (23) will, therefore, give different values, 
and for the two wires. Since the wires are thermal contact the 
tip, their tip temperatures must identical. order establish this 
equality heat will flow from the wire that the better conductor into the 
poorer conductor and the tips will assume temperature that intermediate 
between and the temperature differences being inversely proportional 
the conductivities the wires, thus: 

is, course, consider the use thermocouple with insulation the tip. 


But insulation simplifies testing; and the impression that insulated couples are more 
reliable, although does not wish this point. 


—_ 
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Solving this equation gives 


For purposes illustration constantan may assumed have con- 
ductivity one-twentieth that copper. Consider 
thermocouple which both wires are identical size and insulation, the 
values being Case preceding. Then 0.351° for copper (from 
Table IV), and 1.00. 


For the constantan wire 16,000; 0.01118; 


evident from the foregoing that, when the thermocouple wires differ 
greatly conductivity, the temperature the junction approaches very 
closely that calculated for the better conductor alone. conductivity 
here meant the combined effect specific conductivity and cross sectional 
area. 


Conclusion 


The error indicated thermocouple temperature due the presence 
the insulation may eliminated arranging that the wires follow 
isothermal path for distance from the junction that need not greater 
than the critical distance. When the path the wires not isothermal the 
resulting insulation error determined the variations the ambient 
temperature only over the critical distance from the junction. The magnitude 
this critical distance determined the relative capacities the insulation 
and the wire metal conduct heat. The conditions for minimum insula- 
tion error are small wire diameter, low specific conductivity wire metal, 
thin insulation, and high specific conductivity the insulation material. 


The calculations involved plotting much the curve wire temperature 
are quite laborious, but are, fortunately, not usually necessary. The calcula- 
tion insulation error only relatively simple process, the necessary steps 
being follows: 


(a) Calculate the critical distance Equation (26); 


(b) From the geometry the installation, calculate the ambient tempera- 
ture the critical distance from the junction; 


(c) Calculate insulation error, 


4 
q 
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The above procedure assumes that the ambient temperature curve 
approximately linear over the critical distance. greater accuracy desired 
the procedure would be: 


(a) Calculate the critical distance Equation (26); 


From the geometry the installation, calculate and plot the ambient 
temperature curve for little more than the critical distance from the 
junction; 


(c) With planimeter counting squares, determine the area under the 
ambient temperature curve for the critical distance from the junction; 


(d) Calculate insulation error, T.ds. 


The calculated insulation errors the examples used this paper have 
been relatively small since they were based cases which insulation 
thickness and wire diameter were about ideal and since they were fora 
moderate temperature gradient 50° per inch the tube wall. Reference 
Table will show that gradients times steep may encountered, 
giving under the best conditions insulation errors the order 
degrees. Although insulation error usually not critical the geometric 
errors described the previous paper (1), can excessive adequate 
precautions are not taken. 


hoped that, the future, authors reporting the results heat transfer 
film coefficient measurements will calculate and report the values geometric 
and insulation errors that the accuracy their data will established. 


Table Nomenclature 


diameter thermocouple lead wire. 

term correct for heat entering end area wire. 

specific conductivity wire insulation. 

specific conductivity wire metal. 

distance from center tube any given point the wire. 


distance from surface tube same point the wire defines 
see Fig. 


quantity heat flowing across any given cross section the wire. 
quantity heat flowing into wire through its end area. 


radius the tube. 


distance measured along the thermocouple wire. 


Y 


K = 
M 
1 
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length increment wire; more specifically, “critical distance” 


temperature any given point the thermocouple wire, i.e. 
temperature.” 


average value for any given increment the wire. 


wire temperature the tip junction; choice temperature 
datum also equal the insulation error. 


average temperature the surrounding tube metal which any 
given cross section the wire subjected, i.e., tempera- 


datum temperature. 


ambient temperature distance from the tip 
junction. 


temperature any point the periphery the wire. 


temperature difference across the insulation the wire any point 


its periphery; AT; average values around the periphery 
any given cross section. 


distance measured across the diameter the wire from the hotter 
side. 


thickness wire insulation. 
temperature gradient the tube wall. 
coefficient defined Equation (10a) 


angle defining position point periphery wire, see Fig. 


and indicate different values for two dissimilar metals forming 


thermocouple, e.g., for copper and for constantan. 
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